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Abstract 



A measurement of the differential cross-section for the inclusive production of isolated prompt photons in pp collisions at a center- 
of-mass energy \[s - 1 TeV is presented. The measurement covers the pseudorapidity ranges \rj\ < 1.37 and 1.52 < I77] < 2.37 in 
the transverse energy range 45 < E T < 400 GeV. The results are based on an integrated luminosity of 35 pb~', collected with the 
ATLAS detector at the LHC. The yields of the signal photons are measured using a data-driven technique, based on the observed 
distribution of the hadronic energy in a narrow cone around the photon candidate and the photon selection criteria. The results 
are compared with next-to-leading order perturbative QCD calculations and found to be in good agreement over four orders of 
magnitude in cross-section. 
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The production of prompt photons at hadron colliders pro- 
vides means for testing perturbative QCD predictions [ 1 1, pro- 
viding a colorless probe of the hard scattering process. The 
measurement of the inclusive production of prompt photons 
could be used to constrain the parton distribution functions; 
in particular it is sensitive to the gluon content of the pro- 
ton through the qg — > qy sub-process, which at leading- 
order dominates the inclusive prompt photon cross-section at 
the LHC. 

ATLAS has recently published a measurement of the inclu- 
sive photon cross-section in pp collisions at = 7 TeV using 
an integrated luminosity of 880 nb~ 1 1 3 1 ; a similar measurement 
has been performed by the CMS collaboration |4| using an in- 
tegrated luminosity of 2.9 pb~' . Analogous measurements have 
been perfomed in pp collisions at a lower center of mass at the 
Tevatron [5 6], and in deep inelastic ep scattering at HERA 
[7 8]. This letter presents the measurement of the differential 
production cross-section of isolated prompt photons with trans- 
verse energies Ej above 45 GeV using 34.6 ±1.2 pb~' of pp 
collision data at = 7 TeV collected in 2010. Isolated prompt 
photons in the pseudorapidity ranges \rj\ < 0.6, 0.6 < \rj\ < 1.37, 
1.52 < I77I < 1.81 and 1.81 < \t]\ < 2.37 are studied @. 

In the following, all photons produced in pp collisions and 
not coming from hadron decays are considered as prompt: they 
include both direct photons, which originate from the hard sub- 
process, and fragmentation photons, which are the result of the 
fragmentation of a colored high-px parton ifTUl [TP . Isolated 
photons are considered: from a theoretical perspective, photons 
are isolated if the transverse energy E 1 ™, within a cone of radius 
R = ^{Arf} 2 + (A<p) 2 — 0.4 centered around the photon direc- 
tion in the pseudorapidity (77) and azimuthal angle ((f)) plane |9 1, 
is smaller than E™ 1 . In Jetphox iflOl . used for next-to-leading 



order (NLO) calculations, Ej° is calculated from all partons. 
Similarly, a corresponding isolation prescription is applied ex- 
perimentally on the reconstructed objects, based on the energy 
reconstructed in an R = 0.4 cone around the photon candidate, 
corrected for the effects associated with: the energy of the pho- 
ton candidate itself, the underlying event and the collision pile- 
up 0. The main background to these isolated prompt photons 
is composed of photons from decays of light neutral mesons, 
such as the 7r° or rj. 

Photons are detected in ATLAS by a lead-liquid Argon sam- 
pling electromagnetic calorimeter (ECAL) with an accordion 
geometry, divided into a barrel section covering the pseudora- 
pidity region \t]\ < 1.475 and two end-cap sections covering the 
pseudorapidity regions 1.375 < I77I < 3.2. It consists of three 
longitudinal layers. The first layer has a high granularity along 
the 77 direction (between 0.003 and 0.006 depending on rj, with 
the exception of the regions 1.4 < | ^7 1 < 1.5 and \rj\ > 2.4), 
sufficient to provide an event-by-event discrimination between 
single photon showers and showers coming from a n° decay. 
The second layer has a granularity of 0.025 x 0.025 in 77 x (p. A 
third layer is used to correct for the leakage beyond the electro- 
magnetic calorimeter for high-energy showers, while in front 
of the accordion calorimeter a thin presampler layer, covering 
the pseudorapidity interval I77I < 1.8, is used to correct for the 
energy absorbed before the calorimeter. 

The ECAL energy resolution is parametrized as cr(E)/E = 
a I ~\JE (GeV) © c with the largest contribution coming from the 
sampling term a, corresponding to approximately 10% (20%) 
in the barrel (endcap) region. For energies above 200 GeV the 
global constant term c, estimated to be (1.2 ± 0.6)% ((1.8 + 
0.6)%) in the barrel (endcap) for the 2010 data, starts to domi- 
nate [12|. In front of the electromagnetic calorimeter the inner 
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detector allows the reconstruction of tracks from the primary 
pp collision point and also from secondary vertices, permitting 
an efficient reconstruction of photon conversions in the beam 
pipe and inner detector up to a radius of ~ 80 cm. Further de- 
tails of the inner detector, the electromagnetic calorimeter and 
the whole ATLAS detector are documented in Ref . |[T3l . 

Event samples simulated with Pythia 6.4.21 [14] are used 
to study the characteristics of signal and background events. 
To estimate systematic uncertainties related to the choice of 
the event generator and the parton shower model, alternative 
samples are generated with Herwig 6.5 fl5l . Events used in 
this analysis are triggered using a single-photon trigger with 
a nominal transverse energy threshold of 40 GeV. The trigger 
efficiency, £ tng , is measured using a bootstrap method to be 
(99.4^1)% for prompt photon candidates with E T > 45 GeV 
passing the selection criteria presented below. The same trig- 
ger condition was used for the whole dataset, even though the 
mean number of events per collision rose from < 1 to ~ 3 
as the instantaneous luminosity increased during 2010. Col- 
lision candidates are selected by requiring a primary vertex 
with at least three associated charged particle tracks, consistent 
with the beam interaction region. The total number of selected 
events in data after these requirements is almost 1.7 million, 
with a negligible amount of non-collision background. 

Photon candidates are formed from clusters of energy de- 
posits reconstructed in the electromagnetic calorimeter |16|. 
Clusters without matching tracks are classified as unconverted 
photon candidates. The presence of one or two tracks coming 
from a conversion vertex is used to distinguish converted pho- 
tons from electrons. Converted photon clusters are rebuilt with 
a wider size in <f>, to account for the opening angle between the 
conversion products due to the magnetic field. A specific en- 
ergy calibration IfTBI is then applied separately for converted 
and unconverted photon candidates to account for energy loss 
in front of the ECAL and both lateral and longitudinal leakage. 
Photon clusters are removed if their barycenter lies in the tran- 
sition between the barrel and endcap regions of the electromag- 
netic calorimeter, corresponding to 1.37 < < 1.52, where 
larger uncertainties related to the efficiency measurement are 
expected. Clusters containing cells overlapping with the small 
number of regions with problematic calorimeter readout or with 
very noisy cells are also removed. Over 0.8 million photon can- 
didates with Ej > 45 GeV remain in the data sample. 

A measurement of the transverse isolation energy £™ is as- 
sociated with each photon candidate, computed by summing 
the calorimeter energy in a cone of R = 0.4 around the candi- 
date, as detailed in Ref. [3 |. Corrections to this isolation energy 
are derived from simulation to remove the energy of the pho- 
ton itself that leaks into the isolation cone. An event-by-event 
correction IfTTl [T8l is applied to subtract the estimated contri- 
butions from the underlying event and in-time pileup (i.e. from 
additional proton-proton interactions). The correction to E l ™ is 
typically 900 MeV. After this subtraction, the remaining fluc- 
tuations are dominated by electronic noise from the calorime- 
ter measurement. The effect of the out-of-time pileup, associ- 
ated with collisions taking place in previous bunch-crossings, 
is found to be minimal (i.e. shifts of 200 MeV at most, towards 



lower isolation energies). The corrections mentioned above al- 
low to be directly compared to parton-level theoretical pre- 
dictions. 

All photon candidates having reconstructed isolation energy 
< 3 GeV are considered as experimentally isolated. This def- 
inition is similar to applying a 4 GeV cut on the particle-level 
isolation, defined as the transverse energy of all stable parti- 
cles in a cone of radius R = 0.4 around the photon direction 
(with the underlying event removed as before). The small dif- 
ference between the two, caused by noise and other detector ef- 
fects, is taken into account in the uncertainties associated with 
the photon reconstruction efficiency s rec ° discussed below. The 
particle-level isolation can in turn be related to the parton-level 
isolation in Jetphox that is used for the NLO predictions. The 
efficiency of the isolation criteria is found to be similar (i.e. 
within a few percent) at both the particle-level and the parton- 
level for simulated photons passing the selection described be- 
low. 

As in Ref. Q, the reconstruction and preselection efficiency 
g reco j s com p U ted from simulated prompt photons as a func- 
tion of the true photon Ej. It is defined as the ratio between 
the number of photons reconstructed in a given \rj\ interval 
with reconstructed Ej° < 3 GeV, and the total number of true 
prompt photons with true pseudorapidity in the same \t]\ inter- 
val, and with particle-level transverse isolation energy < 4 GeV. 
The estimated e reco for photons with 45 < E T < 400 GeV is 
~ 85% (75%) in the barrel (endcap) region. The main ineffi- 
ciency (~ 10%) is due to the acceptance loss originating from 
a few inoperative optical links in the calorimeter readout. A 
similar reduction is caused by the isolation requirement in the 
pseudorapidity region 1.52 < I77I < 1.81 where the calorimet- 
ric isolation suffers from larger detector effects. The systematic 
uncertainty on e 1600 associated with the experimental isolation 
requirement is evaluated from the prompt photon simulation by 
varying the value of the isolation criterion by the average differ- 
ence (~ 500 MeV) observed for electrons from W — > ev events 
in data and simulation. The estimated uncertainty varies be- 
tween 3% and 4% depending on 77. The uncertainty associated 
with the imperfect knowledge of the material in front of the 
ECAL is estimated by comparing the expected efficiencies in a 
sample simulated with the nominal ATLAS setup, and one with 
increased material. It varies between 1% and 2.5%, depending 
on 77. 

Shape variables computed from the lateral and longitudinal 
energy profiles of the shower in the calorimeters are used to 
discriminate signal from background lfl6l [191 . As detailed in 
Ref. [3], selection criteria on these variables, optimized in- 
dependently for unconverted and converted photons, are ap- 
plied to reconstructed photon candidates. The requirements on 
these variables are applied in stages resulting in tight candi- 
dates: firstly jets are removed whilst still keeping a high photon 
efficiency and then secondly wide or closely spaced showers 
(i.e. those consistent with jets or meson decays) are rejected. 
The selection criteria have been revised to minimize the sys- 
tematics on the efficiency extraction, especially in the region 
1.81 < \rj\ < 2.37. The photon identification efficiency e 10 
is computed from simulation as a function of transverse en- 
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ergy in each pseudorapidity region. It is defined as the effi- 
ciency for reconstructed (true) prompt photons, with measured 
Ej° < 3 GeV, to pass the identification criteria mentioned 
above. 

Following the same method as Ref. [3|, the value of e ID is 
determined after correcting the simulated shower shapes for the 
observed average differences with respect to data. In the present 
analysis, however, the corrections are estimated for unconverted 
and converted photons separately. This helps to reduce the sys- 
tematic uncertainties associated with the correction procedure. 
The value of e ID varies from 90% to 97%, depending on r\ and 
increasing with Ej. The systematic uncertainty on e ID is also 77 
dependent, ranging from 1.5% to 3%, with contributions from: 
detector simulation; background contamination; (un)converted 
photon misclassification; direct/fragmentation photon fraction; 
the choice of different Monte Carlo generators (MC). These un- 
certainties affect the reconstruction and identification efficien- 
cies in a correlated way, and are treated as such in their combi- 
nation. After applying the isolation criterion and the tight selec- 
tion on the shape variables, almost 173,000 photon candidates 
remain in the data sample. 

As in Ref. 0, a two-dimensional-sideband method is used 
to estimate the background contribution from data and to mea- 
sure the prompt photon signal yield. The two dimensions 
are the transverse isolation energy E™° and the quality of the 
photon, defined by whether or not it passes the shower shape 
identification criteria. On the isolation axis, the signal region 
contains photon candidates with E 1 ™ < 3 GeV, while the 
sideband region contains non-isolated photon candidates with 
Ej° > 5 GeV. On the other axis, the signal photon candidates 
are required to pass the tight identification criteria (tight candi- 
dates). Those failing the tight criteria but passing a background- 
enriching subset of these criteria (non-tight candidates) are con- 
tained in the sideband. A typical distribution of E"° for both 
tight and non-tight data is shown in Fig. [T]for photon candi- 
dates with 45 GeV < Ej < 55 GeV in \q\ < 0.6. The non-tight 
distribution is normalized to the tight one above 5 GeV where 
a only small signal contamination is expected. 

Corrections for the signal contamination in the background 
control regions are computed using prompt photon Monte Carlo 
samples. For the tight isolated signal leaking into the non- 
isolated region, these are as large as 17% at high Ej. Smaller 
leakages of up to 6% are expected for the other two back- 
ground control regions. The purity of isolated prompt pho- 
tons measured with this method increases with Ej from 91% 
at E T = 45 GeV to close to 100% at Ej > 200 GeV. 

The main contributions to the uncertainty on the yields come 
from the fragmentation fraction (< 8%), estimated by conserva- 
tively varying the fraction from to 100% in the signal sample, 
and pileup (5%, with fluctuations up to 8% for 1.52 < | ^7 1 < 
1.81), estimated by increasing the correction to Ej° by 50% 
both in data and simulation. This scaling of the correction min- 
imizes the residual dependency of the isolation on the number 
of primary vertices (i.e. pile-up) in data. The other contribu- 
tions to the uncertainty are: correlated background in the two- 
dimensional-sideband regions (<5% barrel and <10% endcap, 
Ej dependent), definition of the two-dimensional-sideband re- 
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Figure 1 : Distributions of Ej° for photon candidates with 45 GeV < Ej < 
55 GeV in \t]\ < 0.6 passing the tight (solid dots) and non-tight (open triangles) 
shower-shape-based selection criteria. The non-tight distribution is normalized 
to the tight distribution for £™ > 5 GeV (non-isolated region), where the signal 
contamination is fairly small. 



gions (< 5% non-tight and 1% non-isolated), photon energy 
scale (2-8%, 77 dependent), slightly narrower showers in simu- 
lation than in data (2-5%, rj and Ej dependent), isolation shower 
leakage corrections (1-5%), Monte Carlo generator (2%), ma- 
terial effects (< 1%), and prompt electron misidentification 
(~ 0.5%, varying with Ej). Globally, the uncertainties on the 
photon signal yields are less than 10%, and decrease with Ej. 

The average differential cross-section (dtf/dE 1 ™) for the 
production of isolated prompt photons in a bin j of E T ue (in- 
tegrated over one true I77I bin k) is related to the signal yield 
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where e| D '* is the average identification efficiency and /?f is 
the Ej response matrix. The elements of R%. are evaluated from 
the ratio of the true to reconstructed Ej distributions of photon 
candidates, using simulated samples of isolated prompt pho- 
tons. The migration from one Ej bin to another is less than 10% 
in most Ej and 77 regions. A larger migration of up to 18% is 
observed in the region 1.52 < I77I < 1.81, where more material 
is present in front of the electromagnetic calorimeter. Migra- 
tions between rj bins are neglected given the large bin size and 
the excellent ECAL rj resolution. A singular value decomposi- 
tion (SVD) ll20l is used to unfold the Ej distribution for de- 
tector effects. The regularization of the resulting unfolded dis- 
tribution is tuned using simulated events and chosen to be very 
loose to avoid a potential bias toward the truth reference spec- 
trum. The simulation model dependence is tested with pseudo- 
experiments, using Pythia and Herwig simulated samples. The 
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difference of the unfolded cross-section obtained in both cases 
is found to be < 3%. The uncertainty associated with the ECAL 
energy resolution is ~ 1%. The lower and upper Ej constraints 
have negligible effect on the unfolded spectrum. 

The measured inclusive isolated prompt photon production 
cross-sections are shown in Fig. [2] They are presented as a 
function of the photon transverse energy, for each of the four 
considered pseudorapidity intervals. They are also presented in 



tabular form in Appendix A The error bars on the data points 
represent the combination of the statistical and systematic un- 
certainties: systematic uncertainties dominate over the entire 
kinematic range considered. The contribution from the lumi- 
nosity uncertainty (3.4%) is shown separately as it represents a 
possible global change by a common multiplicative factor. The 
data agree with NLO pQCD calculations, obtained with Jet- 
phox 1.2.2 (TO) using the CTEQ 6.6 PDFs EE) and the BFG 
set II [22] fragmentation functions (FF). These predictions are 
negligibly affected when using BFG set I instead. The nominal 
renormalization, factorization and fragmentation scales are set 
to the Ej of the photon. Theoretical calculations using MSTW 
2008 E3 and NNPDF2.0 E) PDFs show a similarly good 
agreement to data. The central values obtained with the MSTW 
2008 (NNPDF2.0) PDFs are 3 to 5% (1 to 4%) higher than those 
predicted using the CTEQ 6.6 PDFs. The total systematic un- 
certainties on the theoretical predictions are represented with a 
solid band. The scale uncertainty (~ 10 %) is the leading theo- 
retical systematic uncertainty. It is estimated from the envelope 
of independent and coherent variations of the three scales, by 
a factor of two around the central value, with the renormaliza- 
tion scale (coherent variation) dominating this envelope at low 
(high) Ej, while the fragmentation scale produces the small- 
est variation. The scale error is summed in quadrature with the 
contributions from the PDF uncertainty (5% at 68% C.L.) and 
the uncertainty associated with the choice of the parton-level 
isolation criterion (2%). The same quantities are also shown in 
the bottom panels after having been normalized to the expected 
NLO pQCD cross-sections. 

In conclusion, the inclusive isolated prompt photon produc- 
tion cross-section in pp collisions at a center-of-mass energy 
sfs = 7 TeV has been measured using 35 pb~' of integrated 
luminosity collected by the ATLAS detector at the LHC. The 
differential cross-section has been measured as a function of 
the prompt photon transverse energy between 45 and 400 GeV, 
in the pseudorapidity ranges 0.0 < I77I < 0.6, 0.6 < I77I < 1.37, 
1.52 < \t]\ < 1.81 and 1.81 < \rj\ < 2.37. In general, good 
agreement between the data and the NLO pQCD predictions is 
observed. This measurement improves the precision and signif- 
icantly extends the kinematic regime explored in the previous 
measurement [3] and is consistent in the region where the two 
measurements overlap. 

Over most of this extended kinematic range the experimental 
errors are smaller than the theoretical ones. The large theoreti- 
cal scale error limits the discrimination between PDFs. Future 
measurements of this process in finer pseudorapidity binning 
and those of the photon + jet system should provide more in- 
sight into the PDF differences. 
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Figure 2: Measured (dots) and expected (shaded area) inclusive prompt photon production cross-sections, and their ratio, as a function of the photon Ej and in the 
range (a) < 0.6, (b) 0.6 < \t]\ < 1.37, (c) 1.52 < \r]\ < 1.81 and (d) 1.81 < |f/| < 2.37. The data error bars combine the statistical and systematic uncertainties, with 
the luminosity uncertainty shown separately (dotted bands). 
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Appendix A. Cross-section measurements 

Tables IA.1BA.4I list the values of the measured isolated 
prompt photon production cross-sections, for the 0.0 < < 
0.6, 0.6 < < 1.37, 1.52 < \rj\ < 1.81 and 1.81 < \rj\ < 2.37 
regions, respectively. The various systematic uncertainties orig- 
inating from the purity measurement, the photon selection and 
identification efficiency and the luminosity are shown. In addi- 
tion, the correlated uncertainties between the efficiency and the 
purity determination are propagated as such and included sep- 
arately (<j corr ). The total uncertainty is the combination of the 
statistical and systematic uncertainties (summed in quadrature), 
except for the uncertainty on the luminosity. 
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Table A. 1 : Measured isolated prompt photon cross-section for \r]\ < 0.6 with statistical and systematic uncertainties. The total uncertainty includes both the statistical 
and all systematic uncertainties (summed in quadrature), except for the uncertainty on the luminosity. 
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iSyield 
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S corY 


^unfolding 
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[pb/GeV] 
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45 


55 


83.3 


0.5 


4.8 


3.3 


3.4 


2.5 


7.2 


2.8 


55 


70 


32.7 
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1.8 


1.2 


1.2 
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1.1 


70 
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12.3 
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0.4 


0.4 


0.4 
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100 


5.3 


0.1 


0.2 


0.2 


0.2 


0.2 


0.4 


0.2 


100 


125 


2.2 


0.05 


0.09 


0.08 


0.07 


0.07 


0.2 


0.07 


125 


150 


0.80 


0.03 


0.03 


0.03 


0.02 


0.03 


0.06 


0.03 


150 


200 


0.26 


0.01 


0.01 


9xl0~ 3 


7xl0~ 3 


8xl0~ 3 


0.02 


9xl0~ 3 


200 


400 


2.8xl0" 2 


2x1 0~ 3 


2x1 <r 3 


lxl0~ 3 


4xl0- 4 


8xl0" 4 


3xl(r 3 


9x10™* 



Table A. 2: Measured isolated prompt photon cross-section for 0.6 < \rj\ < 1.37, uncertainties as in Table|A.l] 
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0.1 
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200 


0.29 
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0.02 


0.01 
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9x1 0~ 3 


0.03 


0.01 


200 


400 


3.2xl0~ 2 
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2xl0~ 3 
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Ixl0~ 3 


4x1 <r 3 


lxl0~ 3 



Table A. 3: Measured isolated prompt photon cross-section for 1.52 < |;;| < 1.81, uncertainties as in Table[A~7] 
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85 
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0.08 
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100 


125 


1.0 
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Table A. 4: Measured isolated prompt photon cross-section for 1.81 < \t]\ < 2.37, uncertainties as in Table|A~Tj 
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100 


125 


1.7 


0.06 


0.1 


0.1 


0.08 


0.05 


0.2 


0.06 


125 


150 


0.55 


0.03 
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